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Abstract 

This paper is a sequel to "Computing diagonal form and Jacobson normal form of a matrix using 
Grobner bases" (Levandovskyy and Schindelar, 2011). We present a new fraction-free algorithm 
for the computation of a diagonal form of a matrix over a certain non-commutative Euclidean 
domain over a computable field with the help of Grobner bases. This algorithm is formulated in 
a general constructive framework of non-commutative Ore localizations of G-algebras (OLGAs) . 
We use the splitting of the computation of a normal form for matrices over Ore localizations into 
the diagonalization and the normalization processes. Both of them can be made fraction- free. 
For a given matrix M over an OLGA R we provide a diagonalization algorithm to compute U, V 
and D with fraction-free entries such that UMV — D holds and D is diagonal. The fraction- 
free approach allows to obtain more information on the associated system of linear functional 
equations and its solutions, than the classical setup of an operator algebra with coefficients 
in rational functions. In particular, one can handle distributional solutions together with, say, 
meromorphic ones. We investigate Ore localizations of common operator algebras over K[x] and 
use them in the unimodularity analysis of transformation matrices U, V. In turn, this allows 
to lift the isomorphism of modules over an OLGA Euclidean domain to a smaller polynomial 
subring of it. We discuss the relation of this lifting with the solutions of the original system 
of equations. Moreover, we prove some new results concerning normal forms of matrices over 
non-simple domains. Our implementation in the computer algebra system Singular:Plural 
follows the fraction-free strategy and shows impressive performance, compared with methods 
which directly use fractions. In particular, we experience moderate swell of coefficients and 
obtain simple transformation matrices. Thus the method we propose is well suited for solving 
nontrivial practical problems. 
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We know what we are, but know not 
what we may be. 

William Shakespeare, "Hamlet". 

1. Introduction 

This paper is a natural extension of our paper (Levandovskyy and Schindelar, 2011). 
We focus on the algorithmic computation of a diagonal form of a given matrix with 
entries in a Euclidean Ore algebra A. In contrast to commutative principal ideal domains 
(Z or for a field K\ there is no such strong normal form as the Smith normal 

form. However, in the case when A is a simple algebra, there is a celebrated Jacobson 
normal form (Jacobson, 1943; Cohn, 1971) of remarkable structure. The existing proofs 
of diagonal and Jacobson normal forms offer an algorithm directly. However, such direct 
algorithms are rarely efficient in general. 

We proposed the splitting of the computation of a normal form (like Jacobson form 
over simple domain) for matrices over Ore localizations into the diagonalization and 
the normalization processes. The diagonalization process can be formulated and applied 
in a fairly general setting, while the normalization depends on the algebra and on the 
properties of concretely given diagonal elements. Both diagonalization and normalization 
can be performed in a fraction-free fashion, and of course, the same ideas apply to the 
computation of the Smith normal form over a commutative Euclidean domain. As we 
demonstrate in Section 6, the question of normal forms over non-simple domains is widely 
open and distinctly differs from the case where Jacobson form exists. 



Email addresses: Viktor .Levandovskyyamath. rwth-aachen.de (Viktor Levandovsk3?y) , 
kristina.schindelarasiemens.com (Kristina Schindelar). 
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In this paper wc; concentrate on the diagonahzation and present the fraction-free 
Algorithm 3.12, which is based on Grobner bases, and its implementation in SlNGU- 
lar:Plural Greuel et al. (2010). As a byproduct, we obtain a fraction-free algorithm 
for the computation of the Smith normal form over a commntativc Euclidean domain. 

Many operator algebras are non-commutative skew polynomial rings (Kredel, 1993; 
Chyzak and Salvy, 1998; Bueso et al, 2003; Levandovskyy, 2005; Chyzak et al., 2007). 
We propose a new class of univariate skew polynomial rings, which arc obtained as Ore 
localizations of G-algebras (OLGAs, see Def. 2.2). This framework is powerful, conve- 
nient and constructive at the same time. Notably, it is possible to perform algorithmic 
computations like Grobner bases for modules in such algebras. Moreover, many compu- 
tations can be done in a fraction-free manner and hence they can be implemented in 
any computer algebra system, which can handle G-algebras or general polynomial Ore 
algebras. 

Such algebras allow to describe time varying systems in systems and control theory 
(Zerz, 2007), (Ilchmann and Mehrmann, 2006), (Ilchmann et al., 1984) and constitute a 
general framework for dealing with linear operator equations with variable coefficients. 
In (Culianez and Quadrat, 2005), applications of a Jacobson form to systems of partial 
differential equations are shown and several concrete examples are introduced. 

In Section 4 we discuss solutions to a given system of operator equations with variable 
coefficients. Starting with a polynomial ring i?* and its 'big' localization R, we denote 
by M the presentation matrix iZ*-module, corresponding to the system of equations. 
We compute matrices U, V and a diagonal matrix D, such that UMV = D . Then, we 
determine a smaller localization R' c i? of i?*, over which both U, V become invertible. 
We provide a module isomorphism for i?'-modules, which allows to treat more general 
solutions, than those one can obtain with the i?-module structure. 

In Prop. 5.1 we propose an algorithm for a probabilistic computation of a Jacobson 
form via a cyclic vector. Moreover, we provide both experimental data from computations 
and analysis of this data. 

At the end, we pose five open problems. In the Appendix, we put the detailed expla- 
nation on the usage of our implementation. 

2. OLGAs and their Properties 

By K we denote a computable field. Describing associative K-algebvas via finite sets of 

generators G and relations R, one usually writes A = K{G \ R) = K{G)/{R). It means 
that A is a factor algebra of the free associative algebra K{G) modulo the two-sided 
ideal, generated by R. Recall the following definition. 

Definition 2.1. Let A be a quotient of the free associative algebra K{xi, . . . ,Xn) by the 
two-sided ideal /, generated by the finite set {xjXi — CijXiXj — dij} for all 1 < i < j < n, 
where c,j e K* and dij are polynomials^ in xi,. . . ,Xn- Then A is called a G -algebra 
(Levandovskyy and Schonemann, 2003; Levandovskyy, 2005), if 

• for all l<i<j<k<n the expression CikCjk • dijXk — Xkdij + Cjk ■ Xjdik — Cij ■ dikXj + 
djkXi — CijCik ■ Xidjk reduces to zero modulo / and 



^ Without loss of generality (of. (Levandovskyy, 2005)) we can assume that dij are given in terms of 
standard monomials . . . x^"' 
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• there exists a monomial ordering -< on K[xi, . . . ,x„], such that for each i < j with 
dij / 0, \m.{dij) -< XiXj . Here, Im stands for the classical notion of leading monomial of 
a polynomial from K[xi, . . . ,Xn]- 

A monomial ordering on K[xi, . . . , Xn] carries over to a G-algebra (and is called ad- 
missible) as in the definition as soon as it satisfies the second condition of the definition. 

As in Levandovskyy and Schindelar (2011), we continue working with Ore localizations 
of G-algebras, which are principal ideal domains. A G-algebra i? is a Noetherian integral 
domain, hence there exists its total two-sided ring effractions Quot(i?) = {R\{0})-'^R, 
which is a division ring (skew field). It is convenient to define a class of subalgebras 
between R and Quot(i?) as follows. 

Definition 2.2. Let i? be a G-algebra, generated by the set X = {xi, . . . ,Xn+i}- 
Suppose, that there is a m with 1 < m < n and a subset Y = {xi^, . . . ,Xi^}, such 
that rfjj.ifc does not involve other variables than those from Y. Then B := K{Y \ 
{xjXi — CijXiXj — dij I i < j,Xi,Xj € Y}) is a G-algebra. For a proper completely 
prime ideal I C B, B/I is a domain. If, in addition, S := {B/I) \ {0}^ is an Ore set in 
A, then we call S~^R an OLGA (Ore-localized G-algebra). It is encoded via the triple 
{R,B,I). 

Prom now on (for simplicity) we consider only OLGAs of the form {R,B,0). 

Prop. 28 of (Garcia Garcia et al., 2009) implies the following characterization of OL- 
GAs. Let Z = X \ Y in the notations as in the above definition. If there exists an 
admissible {Y, Z)-block ordering on R, then B \ {0} is an Ore set in R. 

If we put m = n in Definition 2.2, a corresponding OLGA is a Euclidean domain, 
which we then shorten as OLGAED. Let X \ Y = { d}. In Levandovskyy and Schindelar 
(2011), we proved in Theorem 2.6, that for the case of OLGAED it is enough to require 
the existence of an admissible ordering on A, which satisfies d > Xi^ for all 1 < j < n. 
Moreover, the OLGAED {B \ {0})~^R can be presented as an Ore extension of Quot(S) 
by the variable d. 

2.1. Notations 

In what follows we will work with OLGAED, given as Ore extension R = A[d; a, 5], 
where A = Quot(A*) and A* is a G-algebra in variables {xi, . . . , Xn}. The computations 
will be performed in a G-algebra i?* = A* [ d; a, S] ^ with respect to the monomial mod- 
ule ordering POT (position-over-term), defined as follows. For r, s G Mon(ii*) and the 
canonical basis {efc} of a free module of finite rank, 

re-i < scj i < .j OR (z = j and r < s), 

and r < s with respect to an admissible well-ordering on i?*, in which d is bigger than 
any monomial not involving d. In R, a Grobner basis is computed with respect to the 
induced POT ordering. 



^ Note, that / is an ideal in B and may not be an ideal in R. Moreover, elements of S are identified 
with residue classes modulo /. 

Note, that A, C A C R and A, C C R. 
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3. Praction-free or Polynomial Strategy 

Suppose a matrix M over a non-commutative Euclidean domain R is given. Without 
loss of generality, we suppose that M does not contain a zero row. In this section, we 
show our main approach of this paper. We introduce a method that allows to execute 
Algorithm 3.5 from Levandovskyy and Schindelar (2011) in a completely fraction- free 
framework. The idea comes from commutative algebra (see e. g. (Gianni et al., 1988)). 
Grobner bases were used for the computation of commutative Smith forms in e. g. (Insua, 
2005). 

We define the degree of an element in i?* to be the weighted degree function with 
weight to any generator of A* and weight 1 to d. Thus this weighted degree of / G ii* 
coincides with the degree of f in R and it is invariant under the multiplication by nonzero 
elements in A^,. 

Lemma 3.1. Let M G 7?^^'. Then there exists a diagonal i?-unimodular matrix T G 
A^^^ such that TM G i?*^^"^. Moreover, the computation of such T is algorithmic. 

Proof. If M G i?,^^"^, there is nothing to do. Suppose that M contains elements with 
fractions. At first, we show how to bring two fractional elements ar^b, rr^d for a, c G A^,, 
b,d G R* to a common left denominator, cf. (Apel, 1988). For any /ii, /i2 G A^, such that 
h\a = /i2C, it is easy to see that 

(/iia)-^(/ii6) = a-^h-^^hih = a'^b and (/iia)-^(/i2rf) = (/i2c)-^(/i2rf) = c'^d, 

hence {hia)~^ = a~^h^^ = (/i2c)~^ is a common left denominator. Analogously we can 
compute a common left denominator for any finite set of fractions. Let Tu be a common 

left denominator of all non-zero elements in the i-th row of M, then TM contains no 
fractions. Moreover, T is a diagonal matrix with non-zero entries from A*, hence it is 
i?-unimodular. □ 

Remark 3.2. Note that the computation of compatible factors hi for 01,02 G A^ can 
be achieved by computing syzygies, since {(/ii,/i2) G A^ \ hiai = /12CI2} is precisely the 
module Syz{a\, —02) C A^. The factors hi for more a^'s can be obtained as well. 

Notation. By Q{R,Mt,) we denote the reduced left Grobner basis of the submodule ij.M* 
with respect to the module ordering on i?*, defined in (2.1). Note, that monomials 
of i?*^' are of the form x"'^ ■ ■ -x^^d^ek for ai,/3 G N, 1 < fc < g and is the k-th 
canonical basis vector. Let m be a nonzero vector with entries in i?*. Then by lm(m) we 
denote the leading monomial of m with respect to <* and by Ipos(m) = k the leading 
position of lm(m). By deg(m) we denote deg(lm(m)) = deg(a:;"^ • • • x'^'^d^ek) = P- 

For M G RP^'', we denote by rM = R^^pM the left i?-module, generated by the rows 
of M. 

Define Af* := TM G i?*^' using the notation of Lemma 3.1. Then the relations 
n.,M^ C j^M and /f Af* = rM hold obviously. Thus whenever wc speak about a finitely 
generated submodule rM C R^^'^, wc denote by /jM* a presentation of rM with gener- 
ators contained in i?,. In what follows, we will show how to find i?-unimodular matrices 
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U G i?r^ and V e RV such that 



ri 



U{TM)V ■- 



Since U{TM)V = {UT)MV and UT is a i?-unimodular matrix, our initial aim follows. 



Remark 3.3. Using the fraction-free strategy, two improvements can be observed. On 
the one hand, once we have mapped the matrix we work with from R^^'^ to i?*^"*, the 
comphcatcd arithmetics in tlie skew field of fractions is not used anymore. The other 
improvement lies in the nature of the construction of normal forms for matrices and 
the corresponding transformation matrices. The naive approach would be to apply el- 
ementary operations inclusive division by invertiblcs on the rows and columns, that is, 
operations from the left and from the right. Indeed, there arc methods using differ- 
ent techniques like, for instance, p-adic arguments to calculate the invariant factors of 
the Smith form over Z (Liibeck, 2002), but this method does not help in constructing 
transformation matrices. Surely the swap from left to right has no influence in the com- 
mutative framework. But already in the rational Wcyl algebra _Bi, ^ is an unit in Bi 
and 9^ = — Comparing the multiplication by the inverse element, that is, with 
X, we see that dx — xd + 1 holds. Thus a multiplication of any polynomial containing d 
with the element i in the held of fractions causes an immediate coefficient swell. Since a 
normal form of a matrix is given modulo unimodular operations, the previous example 
illustrates the variations of possible representations. In examples, which we gather in the 
Subsection 7.2, the fraction-free strategy leads to a moderate increase of coefficients. 

On the other hand, switching to the polynomial framework changes the setup. The 
algebra _R* is not a principal ideal domain anymore, which was the essential property for 
the existence of a diagonal form over R. In the sequel, we show how that this problem can 
be resolved by introducing a suitable sorting condition for the chosen module ordering. 
Referring to the argumentation of Lemma 3.5 yields the block-diagonal form 1 with the 
block above. 
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□ 



□ 



(1) 



□ 



Moreover, the rows with the boxed element have the smallest leading monomial with 
respect to the chosen ordering in the corresponding block. A block denotes all elements 
of the same leading position in G{r,M^). In Theorem 3.9 we show that these elements 
indeed generate rM, while in Lemma 3.7 we show that these elements provide us with 
additional information. However, this result requires some preparations. 

Lemma 3.4. Let P he R or R^,. For M e pp^i of full rank and for every 1 < i < q, 
define ai := min{dcg(a) | a G pM \ {0} and Ipos(a) = i\. Then for all 1 < i < g, there 
exists hi G Q[pM) of degree with lpos(/i,) = i. 

Proof. Recall that with respect to d is bigger than any monomial not involving 
d. Let / € pM with lpos(/) = i and deg(/) = ai. Suppose that for all g G Q{pM) 
with leading position i, deg{g) > ai holds. Since Q{pM) is a Grobner basis, there exists 
9 G G{pM) such that \in{g) divides lm(/). This happens if and only if dcg(g) < dcg(/) 
(because i?* is a G-algebra and R is an OLGAED), which yields a contradiction. □ 

The full rank assumption in the lemma guarantees the existence of ai for each com- 
ponent 1 < i < q. Note, that over P = the cardinality of {deg(a) | a G pM\{0} 
and Ipos(a) = i} is greater than one in general, hence there might be different selection 
strategies. We propose to select an element according to min<^ , see Lemma 3.7. Recall 
the Lemma 3.3. from (Levandovskyy and Schindelar, 2011): 

Lemma 3.5. If one orders a reduced Grobner basis in such a way, that \m.{Q{pM)i) < 
■ ■ ■ < lm(^(ijM)m), then [^(i{M)i, . . . ,Q{RM)m] is a lower triangular matrix. 

Corollary 3.6. Lemma 3.4 and Lemma 3.5 yield 

deg{g{RM)i) = min{deg(a) \ ae rM\{0} and Ipos(a) = i}. 
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Lemma 3.7. Let be the degree of the boxed entry with leading position in the i-th 
column, that is 

ai := deg( inin{lm(6) | b € G{r,M^) and lpos(6) = i} ). 
Then for all h G /jM with lpos(/i) =i we have deg(lm(/i)) > ai. 

Proof. Suppose that the claim does not hold and there is /i G rM with lpos(/i) = i of 

degree smaller than ai. By Lemma 3.1, there exists a € A^, such that ah € b„M^,. Then 
deg(a/i) = deg(/i) and lpos(a/i) = i. Due to Lemma 3.4, deg(/) > ai for all / G _r.M* 
with leading position i, hence we obtain a contradiction. □ 

Corollary 3.8. Lemma 3.7 and Corollary 3.6 imply, that for all 1 < i < g 

min{deg(a) | a G rM \ {0} A Ipos(a) = i} = min{deg(a) | a G r.M* \ {0} A Ipos(a) = i}. 

Theorem 3.9. Let M G RP'^p be of full rank. For each 1 < i < p,\et ai be as in Lemma 
3.7. Let us define hi to be the element from {h G G{r,M^) : lpos(6) = i, deg(6) = ai} with 
the smallest leading monomial. Then ji{bi, . . . , bp) = rM. Moreover, the set {bi, . . . , bp} 
corresponds to the subset of all rows with a boxed entry in the block triangular form 1. 

Proof. Since M is of full rank, the minimum in the definition of bi exists for each 1 < 
i < p. Let / G rM\{0}. Due to Corollary 3.8, there exists 1 < k < p such that 
lpos(6fc) = lpos(/) and deg(6A:) < deg(/). Thus there exists an element Sk & R such 
that deg(/ — Skbk) < dcg(6jt). Since / — Skbk G rM, Corollary 3.8 implies that we have 
lpos(/ — Skbk) < lpos(/). Iterating this reduction leads to the remainder zero and thus 
f = J2'!=iSibi. □ 

Notation. Using the notation of the previous theorem, let Q*{jiM) := [bi, . . . ,bg]^ , 
which is by construction a lower triangular matrix. In the sequel, let M G Rp^p be of 
full rank. Then G*{rM) is a square matrix. 

Recall that an involutive anti-automorphism (or an involution) of a ring ^ is a 
ii'-linear map, satisfying 6{ab) = B{h)6{a) for all a,b & A and 9"^ = id^- Moreover, we 
define by d{M) the application of an involution 6 to the entries of the transpose of M. 

Proposition 3.10. Suppose M G Rp^p is a full rank matrix and there is U^, G i?*^^ 
such that U^Mtt = Q{R^Mi,). Let us select the indices 

{^1, . . . C {1, . . . ,4 such that {(f/*M*)t,, . . . , (?7*M*)tJ = G*{rM) (2) 

Then U := [(?7*)t,, . . . , {U^\Y is ii-unimodular in Rp^'p and f/M* = g*(ijM). 

Proof. The equality U M* = Q* (rM) follows by the definition of U. Now we show that U 

is i?-unimodular. Note that f/M, C Rp^'p D M, and r{UM^) = rQ*{rM) ^ rM = rM.^ 
holds. Thus there exists V G Rp^'p such that = V{UM^). Then VU = id^xp and 
analogously UV = idpxp since M has full row rank. □ 

Lemma 3.11. The equality of the following left ideals holds: 

R{e{g*{RM)p^), . . .,e{g*{RM)pp)) = R{g*{e{g*{RM)))pp). 
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Proof. Using the argumentation given in the proof of Lemma 3.4. of (Levandovskyy and 
Schindelar, 2011), wc obtain 

R{o{g*{RMU),...,e{g*{RM)pp)) = R{g{e{g*{RM)))„). 

Bccaiise of R,g*{RM) = RginM) we have 0{g*(RM))R = 9(g(RM))R and thus 
Rg*{0{g*{RM))) = Rg{e{g{RM))). since both gi0{giRM)) and G*{l{g*{RM))) are 
lower triangular matrices, with the latter identity above we obtain R{g{0{g{RM)))pp) = 

R{g*{e{G*{RM)))pp). □ 

Now we are ready to formulate the fraction-free version of the Algorithm 3.5 Diago- 
nalization with Grobner bases from (Levandovskyy and Schindelar, 2011). 

Algorithm 3.12 (Fraction-free diagonalization with Grobner Bases). 

Input: M G RP^P of full rank, 9 an involution on i?* and 9 as above. 

Output: i?-unimodular matrices U,V,D& R^^^ such that U -M-V = D = Diag(ri, . . . , r^). 

Find T e BP^'P unimodular such that TM G ii^^ 

M(o) ^ TM, U ^T, V ^ idpxp 

i ^ 

while M^'^ is not a diagonal matrix or z =2 1 do 

i ^ i + l 

Compute U^"^ so that U^^ ■ M(*-i) = g{R,M('-^)) G i?f ^ 

Select {ii,...,tp} C {1,...,^} as in (2) of Prop. 3.10 

t/«^[(C/«)t,,...,(?7«)t,F 

Af(') ^ 9{G*{rM)) 

if i =2 then 
F^y-^(C/W) 

else 

U ^ • U 

end if 
end while 
return {U,V,M^'^) 

Remark 3.13. It is important to mention, that the matrices U, V, D (hence the elements 
ri as well) have entries from i?*, that is, they are polynomials. However, U and V are 
only unimodular over R and, in general, they need not be unimodular over i?* for obvious 
reasons. In Subsection 7.2 we will investigate, over which subalgebras of R the matrices 
U OT V become unimodular. 

Theorem 3.14. Algorithm 3.12 terminates with the correct result. 

Proof. The proof is a natural generalization of the proof of the Theorem 3.6. from (Levan- 
dovskyy and Schindelar, 2011). At first, we use Proposition 3.10. Moreover, Lemma 3.11 
provides a replacement for the arguments we used in the Lemma 3.4. of (Levandovskyy 
and Schindelar, 2011). □ 

Algorithm 3.12, as well as the original algorithm of (Levandovskyy and Schindelar, 
2011), can be extended to M G i?^'^' along the lines already discussed in Remark 3.7 
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of (Lcvandovskyy and Schindelar, 2011). Our implementation (cf. Section 7) works for 
arbitrary matrices. 

As for examples, a 2 x 2 matrix over the Weyl algebra has been considered in detail in 
Example 3.8 of (Lcvandovskyy and Schindelar, 2011). Note that a fraction-free method 
was used indeed. 

Example 3.15. Consider the first shift algebra i?* = Si = K{x,s \ sx = xs + s) and 
its localization (often called the first rational shift algebra) R = K{x){s \ sx — xs + s) . 
There are precisely two involutions, which can be presented by diagonal matrices on i?* , 
namely x ^ — x, s ^ —s and x i->- —x, s. Let us take the latter and call it 9. Consider 
the matrix in El^^ 



M = 



{x — l)s + x'^ — X xs + x"^ {x + 2)s + a;^ + 2a: 
s -\-x s 



As we can see, T = id2x2 and thus M^'^^ := M, U = id2x2, V = idsxs and i = 0. 
Since M is not a square matrix, under the while condition "while M*^*-* is not a diagonal 
matrix" in the Algorithm we mean the following. The computation will run until the 
matrix we obtain contains a diagonal square submatrix and the entries outside of this 
submatrix are zero. 

1: Since M^^^ is not diagonal, we enter the while loop, i := 1. 



M' := a(fl.MW) = 



-3s2 - (a;2 + 7x + 6)s - - Ax'^ - 3a; {x + l)s^ + (x^ + 2x + l)s 



—3s — 3x 

We set U:=Ui, where f/iM(o) = M' and 



x^ +2x 



Ul : 



s -(a; + 3)s-a;2-4a;-3 
1 -x-2 



Moreover, Af(i) := 0(M') e Rl""^ . 

2: Since M^^^ is not diagonal, we enter the while loop, i := 2. 



-•is^ - x'^ s + bxs + x'^ - Ax'^ + 'ix -3s + 3a; 







-xs — S + X 
x"^ - 2x 




4a;4 + 12a;3 - 4^2 _ 12a; 
— 4a;s — 4s —4a; 

The transformation matrix U2 € Rl^^ is dense, so we show its highest terms with 
respect to s: 



U2 - 



-(a;2 + 5a;- 6)s2 + ... -7(a; - l)s2 + . . . (-a; + l)s2 



-3{x + 6)s2 + . 



-3s^ 



(x2 + 5a; - 6)s2 + . . . 7{x - l)s^ + ... (x - l)s^ 
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Moreover, we put M(2) ;= 6»(M') e -R^x3 and V := e{U2). 
3: Since M^^^ is not diagonal, we enter the while loop, i := 3. 



M' :=e(^,M(2)) = 



+ 12x3 - 4a;2 - 12x -Axs - 4s 
-Ax 



4(x* + 3x^ -x^ - 3x) 
4a; 



[/sM^^^ where [/a 



1 -s 
-1 



Thus we define M^^) — 6I(M') e i?^^^ and U — Us- U. 

4:. Since M'^^) jg diagonal (that is, consists of a diagonal submatrix and the rest of 
entries are zeros) but i = 1 mod 2, we do one more run, which finishes and returns the 
final data: 



D 



[x - l)x{x + l){x + i) 
a; 



U ■■ 



-(x + l)(,x + 3) 

1 -(x + 2) 



-{x'^ + bx- 6)s^ - (x^ +Ax'^ - X- 4)s 

-3(a; + 6)s2 - {x^ + 7a;2 + 5a; - 13)s - - 2a;3 
(a;2 + 5a; - 6)5^ + {2x^ + Gx^ - 14a; + 6)s + a;-* - 



-7{x - 1)8^ - (3a;2 + 5a; - 8)s - 4a; + 4 

- 7a;2 + 8a; - 12 -21s2 _ {Jx'^ + 2x- ll)s - Sa;^ - 2x'^ + 13x - 8 



VIZ 



V23 



- Ix"- 



- 6a; 



7(a; - + (I0a;2 - 16a; + 6)s + Sa;^ - 4a;2 - 5a; + 6 ^33 



where -uis = (x - 1)3^ + {a;^ - a;)s, U23 = Ss^ + (t^ + 2x - 5)s + x'"* - 2x'^ - 3a; + 8 and D33 = 

{-X + 1)5^ - 2(a; - 1)^5 - a;^ + 4a;2 - 5a; + 2 . Indeed, since both nonzero diagonal entries are 
units in K{x)^ the output matrix can be further reduced to 



D' = 



1 
1 



However, then one has to divide explicitly by polynomials in x in transformation matrices. 
Wo arc not going to do this. Moreover, in what follows we will show, how to get important 
information from such matrices containing units from the non-constant ground field. As 
for the concrete example, we conclude, that /B?M = R, thus a system module of M 
is free of rank 1 over R. 



4. Solving Systems of Operator Equations 

Remark 4.1. Let us settle the terminology. Let A be the if-algebra of iiT-lincar op- 
erators. Consider a system of equations in unknown functions wi, . . . ,0Jm- A system is 
linear, if it can be written in the matrix notation, that is 5 • [wi, . . . , Um]"^ = 0, where S 
is a rectangular matrix with entries from A. Then one associates to S' a left j4-moduIe 
M, which is finitely presented by the matrix S. Given a left ^-module T, we usually 
speak of solutions of A4 in J^. The celebrated Lemma of B. Malgrange tells us, that 
the solutions to a linear system of equations S* in a left A-module T arc in one-to-one 
correspondence with the elements of the abelian group Hom^ (7VI , . This allows us to 
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avoid the reference to a specific solution space by adressing an abstract one. However, in 
the context of this paper we address the space of distributions (though yet more general 
hyperfunctions fit into our framework as well) as the space of more general solutions in 
addition to meromorphic functions. 

Assume there is a bigger function space Q, which is an A-module. Thus wc have an 
exact sequence — >■ — ?■ 5. Due to the left exactness of the Horn functor, we obtain 
that }iomA{M,J^) — > Hom^(A1,^) is an exact sequence as well. This justifies the 
fact, that working with a more general solution space Q we obtain not less solutions to 
a module Ai as with J^. 

An _R-modulc is naturally an _R*-module, but an i?*-module is not necessarily an R- 
module. There can be i?*-modules M with 5-torsion, that is those for which S~^M = 
holds. Hence, though working over R brings significant comfort, we are interested 
in gaining more information from the algebraic structure by looking at i?* and, more 
generally, subalgebras between i?* and R. Since there are i?,-modules, which are not 
i?-modules, in view of Malgrange's Lemma there might be more solutions in i?*-modules 
as in i?-modules. Sec Examples 4.9 and 4.10 for illustrative details. 

We need to recall and develop mechanisms, which allow us to tackle localizations 
of operator algebras. We will show, how different small enough localizations look for 
common operator algebras. 

4-1- Ore Multiplicative Closure 

Let W e R'^^'^ be a square matrix with entries in R. Assume that there exists left 
inverse matrix T, that is TW = idr-xr- Then by Lemma 3.1 there exists a diagonal 
matrix Q = Diag(. . . ,qii, . . .) such that Q (resp. QT) has entries from A^, (resp. i?*). Let 
^ = {qu I 1 < i < (Zii ^ K} U {1}. Denote by Rw the Ore localization of i?* with 
respect to a multiplicatively closed Ore set Sw, which is defined as follows. Let A4(f2) be 
the two-sided multiplicative closure of ft, equivalently the (possibly non-commutative) 
monoid; generated by a finite set f7 C i?* \ {0}. Let Sw be an Ore closure of A4{n) in 
i?*, that is a set, containing M{il), which is an Ore set in Such a closure always 
exists since i? is a localization of i?* with respect to 5 = ^4* \ {0}, but we are interested 
in computing a closure, which is minimal in the sense that there exists no closure S' 
satisfying M{n) C S' C Sw 

4-2. Ore Closure in Classical Algebras 

In order to localize a domain R with respect to a multiplicatively closed set 5, the 
latter needs to be a (left and right) Ore set in R. Thus, wc are going to show, how 
to compute an Ore closure of S in R. Moreover, we ask for small generating sets of a 
monoidal localization. 

Let us recall the well-known Lemma (see e. g. Zariski and Samuel (1975)) first. 

Lemma 4.2. Let / e K[x] := K[xi, . . . be a non-constant monic polynomial and 

^ = {P I * € No}. Moreover, let / = /f ^ • • • ff'' be an irreducible factorization in K[x] 
and T := {f^\. . . , J"- | a G N^}. Then S-^K[x] = T-'^K[x]. 
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Thus, if / = (a; + ify^, then S={{x+ if^y^^ | i G No} and T={{x + l)\ y^\i& No}. 

Let 5 be a multiplicatively closed subset of a domain R and ^ S. In order to prove, 

that S is an Ore set in i?, one has to show, that V (s, r) E S xR, there exist {t, q) E S x R, 
such that r ■ t = s ■ q in R and V(i, q) £ S x R there exist (s, r) & S x R, satisfying the 
same condition. In other words, one can rewrite any left (resp. right) fraction as a right 
(resp. left) fraction. 

Wc will analyze the smallest nontrivial Ore sets in the first Weyl, shift and q-commutative 
algebras (these results seem to be folklore) and give simple constructive proofs of the Ore 
property for them. 

Lemma 4.3. Let Ai be the first Weyl algebra and / G K[x]\K. Then S = {f \ iG No} 
is an Ore set in Ai. 

Proof. We have for and i G N, that d • = • (/d + (z + l)-g^). By induction, 
one can prove, that for j G N and i + 1 > j one has d-' • = ■ (/■? d^ + vu), 

where the terms of Vij G Ai have degree at most j — 1 and contain derivatives up to f^^\ 
Suppose we are given g = Ylj^o bj {x) G Ai with ^ and f'' for a fixed fc G N. 
Then 

Thus 

d d 

□ 

Consider now the first shift algebra 5*1 (cf. Example 3.15). Since sx = {x + l)s, we see 
that for all ^ G Z 

{x + z + i)~^s = s{x + 

thus it seems natural to have polynomials with integer shifts of their argument in the set 
S as above in addition to S itself. 

Remark 4.4. For / G K{x] \ K, the set 5 = {/* | i G N} is not an Ore set in the 
first shift algebra. Take s and f''{x) G S, we're looking for f^{x) and t G Si, such 
that sf{x) = ,f{x)t. The left hand side is f{x + l^s thus ,f{x)t = f{x + ifs. But 
f{^) \ fi^ + 1) fo'^ non-constant /, thus there exists no < G 5*1 satisfying the latter 
identity. It means, that we have to enlarge S in order to obtain an Ore set in Si. 

Lemma 4.5. Let Si be the first shift algebra (cf. Example 3.15) and / G K[x]\K. Then 
S = {/"(a; ±z) \n,z € No} is an Ore set in ^i. 

Proof. Given g = J2j=o^ji^)^^ ^ '^i ^^^^ ba and h{x) = f^{x + zq) G S with 
G N, ^0 G Z, let us define gf{x) := 11^=0 - i) G 5. Then 



g-gfix)=J2bjix)s' ■U^ix-i) = h{xf-^bjix){ [] hix+j-i)s^). 

j=0 i=0 j=0 i=0,i^j 



□ 
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A similar phenomenon can be observed in quantum algebras as well. 

Lemma 4.6. Let Qi be the first g-commutative algebra K{q){x,s \ yx = qxy) and 
/ e K[x\ \ K. Then S = {/"(9=^^a;) \n,zG No} is an Ore set in Qi. 

Proof. Note, that for any g{x) £ K[x] one has y"^g{x) = g{q"^'x)y. Suppose wc are given 

g = J2j=o ^ji^)y^ S Qi ^■iid H^) = f'^il^x) e Si. Let us define gf{x) := HiLo h{q~^x) e 
S. Then 

d d d d 

j=0 i=0 j=0 i=o,i^j 

□ 

Let ii' be a ii'-algebra with C R' C R, where R^,R are as above and M e -R™^". 
Consider a system of equations Mlo = in unknown functions io — (wi, . . . ,a;,„) from a 
space of functions J^, which will possess some module structure, see below. Assume, that 
we have computed U,V (unimodular over R) and D = Dia,g{dii, ... ,dmm) satisfying 
UMV = D. 

4.3. T is an R-module 

Inverting V, we obtain UM = DV~^ , hence Mco = is equivalent to = UMw = 
DV~^uj. Thus, introducing an i?- automorphism of defined by zu := V~^uj, wc obtain 
a decoupled system {duWi = 0}. Note that if da = 0, then zui is called a free variable of 
the system, e. g. in (Zerz, 2006). The solutions of the decoupled system in are precisely 
the solutions of Mco = in .F. 

4.4- is an Ri, -module 

Indeed, the i?-automorphism of above can be defined as soon as V is invertible. We 
can regard this as a kind of "analytic" transformation of J^. 

Proposition 4.7. Let UMV = D as before. Let Sjj and Sy be Ore multiplicative 
closures of U and V respectively, according to Sect. 4.1. Moreover, let S be an Ore 
multiplicative closure of the monoid Sy U Su- Then, over S'^^i?* C i? we have Mw = 
<S4> D{V~^uj) = 0. Thus it is possible to decouple the system M. Moreover, there might 
be solutions in an S'~^i?*-module G- 

Proof. Let Rv = {Sv)~^R* ■ Then on an i?v-module Q we can define an i?v-automorphism 
w V~^io. Thus over Ry we have UMui = <^=^ Dw = 0. Moreover, U is invertible over 
{Su)~^R*. Hence, over S~^R» we have Moj = <^ UMoj = <^ Dw = 0. □ 

Corollary 4.8. With notations of the Proposition, there is an explicit isomorphism of 
5" ^i?- modules 

(5-^i?)"/(5-^i?)"M ^ (5-^i?)™/(5-^i?)"C/My = 

Note that the left transformation matrix U also can contain essential information 
about the so-called singularities of a system, which often are connected to meromorphic 
(and non-holomorphic) solutions. 
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Example 4.9. Over the first Wcyl algebra Ai = i?*, consider ttic single equation 
{xd)oj = 0. Here U = x and thus Su = {x' \ i £ No}. So Ru = {Suy^R* = 
K[x,x~^]{d) C K{x){d) and over Ru, the equation is equivalent to doj = 0, whose 
solutions in any nonzero i?-module contain K. 

On the other hand, the division by x over i?* is not allowed. Consider an i?*-mod\ile 
I?'(R) of distributions. Then, we see that a if- multiple of the Heaviside step function is 
a solution. Thus, K ■ H{x) ® K \s & subspace of solutions of the system (a; d)w = 0. 

Example 4.10. Consider the univariate sequence space 5, that is the vector space 
of all functions f : Z ^ K, which is an S'i(iir)-module. Recall, that a discrete analogue 
H{n) of the Heaviside step function is defined to be 0, if n < and 1 otherwise. 

Consider the analogon of the equation above (n(s — l))w = over the shift algebra 
Si{K) = J?*, cf. Example 3.15, where s acts as {suj){n) = u){n + 1). The invertibility of 
n is reflected in the localization with respect to {{n ± fc)™ | fc, m 6 No} (cf. Lemma 4.5) 
and implies that there is 1-dimensional subspace of constant solutions. 

Over J?*, n is not invertible. Moreover, the ideal (n) C i?* annihilates any constant 
multiple of the Kronecker delta (5„^o £ S. Since (5,i_o = H{n) — H{n—1) = (s— l)(i7(n— 1)), 
another set of solutions to the considered equation are if-multiples of H{n — 1). Hence, 
as one can easily check, K ■ H{n — 1) ® K c S is indeed the whole set of solutions to the 
equation (n(s — l))a; = 0. 

5. Cyclic Vector Method 

The existence of the Jacobson form of a matrix over a simple Euclidean domain (Cohn, 
1971; Jacobson, 1943) is a very strong result. In particular, for a square matrix M of full 
rank over R, a Jacobson form is Diag(l, . . . , 1, r) for some r £ R \ {0}. Then a module, 
presented by M is isomorphic to a cyclic module and its presentation is a principal ideal. 
The method of finding a cyclic vector in a finitely presented module and obtaining a left 
annihilating ideal for it is also used in D-module theory. 

Proposition 5.1. Let i? be a simple OLGAED, representable as A[d: a, S] with a division 
ring A over the field K. Let M = Diag(mi, . . . , rrir) be a full rank r x r matrix, that is 
rrii ^ 0. Then d := d{M) = ^deg(mi) is an invariant of the module R^/R^M, since it 
is the dimension of the module over the division ring A. Let p = [pi, . . . 6 R^^^ 

and c G i? be a generator of the left annihilator ideal of p in the module R^/R^M. If 
deg(c) = deg(mj), then Diag(l, . . . , 1, c) is a Jacobson form of M. 

Proof. There is an i?-module homomorphism ipp and the corresponding induced exact 
sequence 

— )■ R/{c) = R/kenpp R'/R^'M — y cokercpp — y 0. 

Since all i?- modules above are finite dimensional over A, the dimension of cokerc^p is 
precisely d — deg(c). Hence if deg(c) = d, then 

i?7i?'' Diag(l, . . . , 1, c) ^ R/{c) ^ R''/R''M = R^'/R' Diag(mi, . . . , m^). 

□ 



15 



Remark 5.2. By using the previous Proposition, we propose the foUowing probabilistic 
approach for the computation of a cychc presentation of a module. We use the dimension 
d of M as the certificate. For every 1 < i < r, consider a polynomial pi of degree at most 
deg(mi) — 1 in d with random coefficients from A. Compute the generator c £ i? of ker (pp. 
If deg c = d,we are done. Otherwise (it means that the image of (pp is a proper submodule) 
one takes another set of random polynomials pi and repeats the procedure. In order to 
turn this approach into an algorithm, one needs to obtain probabilistic estimations on 
the length of random coefRcients like in (Kaltofen et al., 1989) and (Storjohann and 
Labahn, 1997). However, it is more complicated in the case when A is non-commutative 
division ring and needs to be investigated in more detail. Our experiments, illustrated 
by the examples below, detect a considerable coefficient swell both in the intermediate 
computations and in the output. Thus we conclude, that the deterministic method by 
Leykin (Leykin, 2003) is better for the case, when random vectors contain polynomials 
from K[x] as coefficients. 

Yet another application of this approach is the search for a proper submodule of a 
module. Namely, if deg c< d, the image of ipp is such a submodule. 

In the following examples we work in the first rational Weyl algebra A. Consider 
two 3x3 diagonal matrices with polynomials of degrees 1,2,3 in d at the diagonal. 
We call these matrices Mi and M2. By the remark above, in order to find a cyclic 
vector of the corresponging module, it is enough to consider a vector p of the form 
[c{x), a{x) d + b{x), u{x) d? + v{x) d + w{x)^ for a, 6, c, u,v,w € K[x]. We generate the 
latter polynomials in a random way, taking their coefficients from Q from the range 
[0, . . . , 100). The vector pi has coefiicients of a;-degree at most 3 and the vector p2 of 
a;-degree at most 4: 

Pl = [98x^ + 4, {2x^ + 17) d + 87x^, (98x^ + llx) + {8x^ + 62x^ + 31) d + 89a;]'^ , 

P2 = [50x4 ^ j^3^3 _,_ 97^2 _ (-25^.4 g-^^ + 72) d + SSx-* + 96x^ + 90x, {41x^ + 57x^)(f + (36x^ + 53x^ + 

54a; + 83)d + 2a;3 + 87]^. 

By computing kernels as in Prop. 5.1 we obtain generators c^, C2 for i = 1, 2 respec- 
tively. Since degc* — 6 — '^^^i deg Mkk, it follows that pj are cyclic vectors for both 
matrices Mi, M2 and the Jacobson form of Mj with respect to pj is Diag(l, 1, c*). 

Example 5.3. Consider the matrix Mi = Diag(d,a;d^ + 2d, x^d^ + 4a:; d^ + 2d). For 
the vector pi we obtain the cyclic generator c\ = (1011752a:* - 348435a:^ - 846320a:5 - 
2965480a;'') d« + (9105768a;^ - 3484350a:« - 101558401* - 38551240a:3) d^ + (15176280a;8 - 6271830a:5 - 
25389600.r3 - 115653720a:2) d* - 35585760a: d^ + 35585760 d^. 

For the vector p2 

we obtain the cyclic generator = (395647200a;i6 _ 264348000a;i5 + 72499244880a:i'* - 

432343555560x1=* + 267049574380x12 + 8987298499008x" - 45076322620512x-i'' + 91959270441432x^ - 
24315286945590x« + 20084358713472x^ - 19034034270714x'5 - 3963517931016x5 + 240924562515x"' + 
239689051938x3 )d6 + ... . 

Example 5.4. Consider the matrix M2 = Diag( d, a; d^ + 2 d, x"^ d^ -|- 3a; d^ -|- d). For the 
vector Pl we obtain the cychc generator cf = (8352x" +149292x11 +32l3954xi° -8623701x8 + 
46759968X* + 18251056x'^ + 12525848x8 + 1085484x5 + 968320x4) d<^ + . . . . 

For the vector p2 we obtain the cyclic generator c| = (544946x^1 -8586000x20+4018843296x1^ - 

2355826816422x18 +7580096636636x" -91394273346228x1*' +394488979119486x15 - 1039358677414560x14 + 
2049350822715951x13 -6702303668155704x12 +15453420608607570x11 -9398963461913820x1° -1453404219438726X''- 
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83615730159492X* +313903596685743I;'' +646008054793470x8 +64626268386222a;B -13528366907400a;*- 
3553905268080x3 )d6 + . . . . 

Analysis of the data. Since the computed cycHc generators are fraction- free, we 
present them as polynomials from Z[a:, d]. We determine the following data and put 
them into the table. 

• NT number of terms, TD total degree; 

• BC, SC and AC stand for the biggest rcsp. the smallest resp. the arithmetic average 
of absolute values of coefficients; 

• BX, SX and AX stand for the biggest resp. the smallest resp. the arithmetic average 
of degrees in x in every monomial a;" d*". 



Poly 


NT 


TD 


BC 


SC 


AC 


BX 


SX 


AX 


^1 


14 


14 


1.15-10^ 


3.48-105 


2.14-10^ 


8 





4.3 




85 


22 


3.2-1015 


2.6-10^ 


2.3-10^^ 


16 





6.8 




43 


18 


3.9-10^ 


8.3-103 


6-10^ 


12 





5.8 




126 


27 


4.95-10" 


5.4-10^ 


2.35-10^6 


21 





9.4 



Let us see, what can happen when the random polynomial coefficients from above are 
taken to be just numbers from K. 

Example 5.5. For the matrix M2, let us take 5- and 10-digit random integers. The 
corresponding vectors are then pi = [5535, 3892 d + 20690, 6069 + 20660d + 17323]^ and 
P2 = [1109725034, 618308146 d + 1684065511, 2034108815 + 1702526110d + 1361184996]T . The 
corresponding cyclic generators arc then ci = (5l973949l8lla;i'^ + . . .) d"* + . . . and C2 = 
(689281531286620056404507706110 + . . .) d"' + . . . . 

The coefficient data is put into the table with the same notations as before. 



P()l\- 


NT 


TD 


BC 


SC 


AC 


BX 


SX 


AX 




20 


14 


2.14-1()J-' 


4.77-l()i" 


.■).48-l()i- 


10 


:-! 


(),5 


C2 


20 


14 


6.15-10^^ 


6.89-10^6 


8.29-10^^ 


10 


3 


6.5 



As we can see from the degree in d. both ci, C2 are not cyclic vectors. Notably, numer- 
ous experiments in the setup of this example never led us to a cyclic vector for A'' / A''M2. 
Instead we repeatedly obtained polynomials Ci of degree 4. This shows, that there is a 
4-dimensional submodule N of the 6-dimensional module /j^M2-, which dominates 
over others for the special choice of the form of test vectors. 

6. Normal Form over a General Domain 

As we have noted, there are many different-looking normal forms over a non-simple 
domain. 
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Definition 6.1. Let i? be a ring. An element r E R is called two-sided, if r is not a 
divisor of zero and /j(r) = (r)/{. It is called proper, if ji{r) C R. 



Though in Cohn (1971) a two-sided element has been called invariant, we propose 
to use two-sided instead, due to the ubiquity of the word invariant. 

Let i? be a domain and a /S'-algebra. Then r G R is proper two-sided if and only if 
V.s, s' E R3t,t' G R such that rs = tr and s'r = rt' . If R admits a Grobner basis theory, 
this is the same as to say "{r} is a two-sided Grobner basis of the two-sided ideal (r)". 
It is straightforward, that for a simple domain is the only proper two-sided element. 

Generalizing the statement in the Example 4.4 of (Levandovskyy and Schindelar, 2011) 
leads to the following result. 

Lemma 6.2. Let i? be a non-simple Euclidean domain and a if-algebra. Moreover, let 
m < n be natural numbers and r G i? be a proper two-sided element. Then for the 

2x2 matrices Di — Diag(r'",r") and D2 = Diag(l, r'"+") the corresponding modules 
Ml = R^/R^Di and M2 = R^/R^D2 are not isomorphic. 

Proof. At first, we note that the set of proper two-sided elements is multiplicativcly 
closed, thus r''\ r"\ r'"'^"'' arc proper two-sided. As we can see, Ann/j Af2 = (r'"+") and 
Anui^ Ml = Annfl(i?/(r'")®i?/(r")) = (r'")n(r") = (r"). We have to show that the two- 
sided ideals Ann/j Mi, Annji M2 are not equal, then the claim follows. Let S = R®kR°^^ 
be the enveloping algebra of R and l : R ^ S & natural embedding of iiT-algebras, 
then any two-sided ideal of i? is a left ideal over S. Suppose that for two-sided ideals 
(r") C (r'"+"). Over S we have the inequality of left ideals (t(r")) C (t(r™+")), hence 
there exists s £ S*, such that ;,(r") = S6(r™+") = si(r")t(r"). Then (st(r™)-l)-t(r") = 
and hence st(r'") = 1, since S* is a domain. Thus, (t(r'")) = S* as a left ideal and 
hence (r™) = R, what is a contradiction to the assumption that r is a proper two-sided 
element. □ 

Recall that for a, 6 £ i?, one says that a totally divides h (and denotes a || 6) if and 
only if there exists a two-sided c € R, such that a | c | 6^ (Cohn, 1971). Moreover, a\ \ a 
if and only if a is two-sided element. 

Example 6.3. In the first shift algebra ii* = Si, cf. Example 3.15, consider the diagonal 

matrices 



and D2 



s-\-x 
s{s + x) 



s{s + x) 



At first, let us analyze the appearing elements for the proper two-sidedness. By computing 

two-sided Grobner bases over 5i, we obtain that s is already such a basis, whereas 
Si (s + x)si = Si {x, s)si and Si {s{s + x))si = Si {{x + l)s, s^)si ■ Thus neither s + x uoi 
s{s + x) are proper two-sided in ii*. In the localization R = {K[x] \ {0})~^i?*, we even 
see that r{s + x)r = R and ji{s{s + x))ji = r{s)r. 

In the matrix D2 there is a total divisibility on the diagonal. Indeed, s is proper 
two-sided and s \ s{s + x). Hence, for any f G R the result of reduction of s with 
fs{s + x) = f{s-\-x + l)s will be (1 — f{s + x + l))s € (s). Hence, no reduction procedure 
will lead to a unit instead of s. 



^ Recall, that a | c if 3d € R such that either c = ad or c = da holds. 
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On the other hand, we see that s{s + x) € (s). Then s{s + x) — fs = {s + x+l — f)s 
has degree 1 in ,s only for f = s — g{x). where g{x) ^ — (x + 1). This reduction produces 
[0, s(s + x)]^ — {s ~ g{x)) ■ [s, sY = [— (.s — g{x))s, {g{x) + x + l)s]"^, which does not 
lower the degree in the column but only interchanges terms of degrccis 1 and 2. Hence 
no further essential simplification is possible, since replacing .s(s + a;) by a similar factor 
can neither lower the degree nor lead to a two-sided element. 

In the matrix Di, there is no total divisibility on the diagonal, hence we will be able 
to achieve the lower degree. And indeed, for 



U = 



1 1 
+ XS + 2S +XS + 2S + X + 2 
which are i?-unimodular, we obtain that 



UDiV 



V = 



s + 1 —s^ — xs — s 
-1 s + x-1 



X 




1 


X + 2 




s(s + x)(s + a; — 1 



Thus Di possesses even the Jacobson normal form Diag(l, s{s + x){s + x— 1)) over R. 



7. Implementation and Examples 

Our implementation of the computation of a diagonal form together with transforma- 
tion matrices is called jacobson. lib. It is distributed together with Singular (Decker 
et al., 2011) since its version 3-1-0. In the Appendix we put an example of a SINGULAR 
session with input, output and explanations. 

7.1. Comparison 

There are other packages to compute diagonal and Jacobson forms. The package 

Janet for Maple (Blinkov et al., 2003; Chyzak et al., 2007) directly follows the classical 
algorithm with no special optimizations. In Maple packages by H. Cheng et al. (Becker- 
mann et al., 2006; Cheng and Labahn, 2007; Davies et al., 2008) modular (Modreduce) 
and fraction-free (FFreduce) versions of an order basis of a polynomial matrix M from 
an Ore algebra A are implemented. The computation of the left nullspace of M and 
indirectly the Popov form of M uses order bases. There are also experimental implemen- 
tations, mentioned in (Culianez and Quadrat, 2005) and (Middeke, 2008). 

In (Levandovskyy and Schindelar, 2011), we compared our implementation with the 
one by D. Robertz. In turned out that with our implementation, one experiences moderate 
swell of coefficients and obtains uncomplicated transformation matrices. 

Unlike in theoretical considerations, our implementation usually returns diagonal ma- 
trices with elements of descending degrees on the main diagonal. 

7.2. Examples 

Example 7.1. Consider two versions of the Example 3.8 in (Levandovskyy and Schin- 
delar, 2011). 
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(A). In the first Weyl algebra = K[x][d;id, consider tlie matrix 

^2 - 1 a + 1 
+ 1 d-x 



M = 



2x2 



Then the algorithm computes 



= A 



-(x + 1)0 + + 2; _|_ 1 (x + 1)9 + 3; 



V = 



(x + 1)92 + 2d-x+l 1 



a - x -a - 1 

Let lis analyze the transformation inatriecis for i?*-uniniodiilarity. Indeed, V is uni- 
modular over i?* since it admits an inverse V. However, U is unimodular over S~^R^ 
for S = {{x + iy \ i G No} (cf. Lemma 4.3) since U ■ Z = W and W is invertible in the 
ring containing the inverse of {x + 1)^, that is over a ring, containing S~^Ri,. 



V 



-(x + l)d? - 2A + X- 




1 1 



2(d + l) {x + l)A + x-2 
2{A-x) {x + l)d-x'^ - x- 



-4(x + 1)2 
2 h{x + 1) 



(B). In the first shift algebra S\, cf. Example 3.15, one has 

Diag( (a; + \){x + 2)s^ + 2(a; + l)s - {x - \){x + 2), 1) 



-{x + l)s + x{x + 2) (a; + l)s + a; + 2 

-s + x + 1 s + 1 



1 s + 1 

\ s — X 



1 



-(x + 2)s2 





2s + X 1 



Denote the equality above as D = UMV ^ it turns out, that V is even i?«-unimodular. U 
becomes invertible in a localization, containing (x+2)~^, thus, by Lemma 4.5, containing 
a ring S'^R^ for S = {{x± n)™ | n, m G No}. 

Example 7.2. Consider now the same matrix as in Example 3.15 over the first Weyl 
algebra -R* = K{x, d | dx = xd + 1), respectively R = K{x){A \ da; = xd+ 1): 



M = 



{x-l)d + x^ ■ 
d + x 



-X xd + x^ {x + 2)d + x^ +2x 



We obtain 

3x^{x + 2)^ 



D = 











,U = 



x{x + 2) d - 2(1 + 1) -x(x + 2)2 d - {x"^ + 4x^ + 3x'^ - Ax - i) 
(2: + 1) d - 2 -{x + l){x + 2)d-{x^ + 3x^ + x- 3) 



V 



-{x* -9x'^)d^ - {x^-x^ ■ 



I8x) d - {6x^ - 24x2 ^ ;^g) 



-{3x^ -27x)d'^ - {x^ + 5x* 



Vi2 

■ 9x^ - 18a;2 - 81) d - (x^ + 2x^ - 61" + Sa:^ + 9a:2 _ 543.) 
(a;" - 9a;2) d2 + {2x^ - lOx^ - 18a;) d + (a;^ - 15a;2 + 18) V23 1 

where Vi2 = {x^ + 3x-2 - 2x-) d2 + (x-^ + Sx-"* + 4x2 + 6x + 2) d + (5^3 _,_ ;^2x2 - 6), H22 = (3x2 + 
9x - 6) • d2 + (x* + 8x3 _,_ 2^33,2 _,_ _,_ 27) . d + (x^ + Sx"* + 6x^ + 14x2 + 42x + 23) and 1^23 = 
-(x3 + 3x2 _ 2x) d2 - (2x4 _,_ g^3 _|_ 2x2 + 6x + 2) d - (x^ + 3x'^ + Sx^ + 18x2 + 12x - 6) . 
U is unimodular in the localization with respect to Su, which is generated as monoid by 
a;,a; + l,a; + 2. l^is unimodular in the localization with respect to Sy, which is generated 
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by x, X + 2, X + 3, X — 3. Thus, the isomorphism of modules lifts from Rto S ^R*, where 
S' = {x*i(x + l)*=(x + 2)*3(a. + 3)u(a._3)«5 eno}. 

The same system over the shift algebra was computed in detail in Example 3.15. There 
it turns out, that U is unimodular over a ring containing (x + 2)~^ and V is unimodular 
over a ring containing the inverses of {x — 2, x — 1, x, x + 3, x + 4}. Due to Lemma 4.5, 
the needed Ore set for the localization is not smaller than S = Su = Sy = {{x ± n)'^ \ 
n,m € No}. 

8. Conclusion, Further ResecU-ch and Open Problems 

8.1. Further Research 

We do not perform an analysis of the theoretical complexity of our algorithm, since it 
has been designed by using Grobner bases, thus the formal complexity is too high. Indeed, 
one can see, that we use Grobner bases rather for convenience, that is in order to present 
an algorithm, working over arbitrary OLGAED. The original (non- fraction- free) algo- 
rithm can be reformulated in terms of cixtendcid left and right greatest common divisors 
over a concrete algebra. To the best of our knowledge, the theoretical complexity of such 
GCD algorithms depends on the given algebra. For rational Weyl algebras over general 
fields of characteristic zero there is a paper (Grigor'ev, 1990). It is interesting to find es- 
timations for classical operator algebras like the shift algebra, q-Wey\ and g-commutative 
algebras with rational coefficients. Having this information, one could perform complex- 
ity analysis for the kind of algorithms we propose. It was reported in (Middeke, 2008), 
that the Jacobson form can be computed in polynomial time. In our opinion this should 
be formulated more precisely in the framework, proposed in (Grigor'ev, 1990). We do 
not claim that our algorithm is superior in terms of theoretical complexity to the others 
but stress, that its fraction-free version is widely applicable in practical computations. In 
particular, our algorithm returns transformation matrices with reasonable sizes of their 
entries and their coefficients. 

Concerning the further development of our implementation in Singular (Decker 
et al., 2011) called jacobson. lib, we plan the following enhancements. In order to pro- 
vide diagonal form computations over (/-algebras like (/-commutative, g-shift, (/-difference 
and (/-Weyl algebras, we need to find involutive antiautomorphisms. It turns out, that 
there are no K (q)-\iiiea,T but if-linear involutions, which act on (/ as a non-identical in- 
volutive automorphism of K{q). We will work on finding such invohitions and adapt the 
implementation to the more general situation. Moreover, it is possible to employ modular 
Grobner bases in the implementation. 

We have demonstrated, how working with a ring ii* and its localization S~^Rt, with 
respect to some Ore set 5 C i?* is used in various situations. In this paper we concen- 
trated on rings of polynomials respectively rational functions. The same ideas in theory 
immediately apply to operator algebras i?* = i4r[[x]](d) and R = if((x))(d). However, 
in practical computations with the computer we are restricted to finite power series, 
presented by a Laurent polynomial. 

For the Weyl and g-Weyl algebras over a field K with chaxK = 0, algorithmic compu- 
tations are possible over the local ring i4r[x]mp for p G Jf" and rrip = (xi — pi, . . . , x„— 
Notably, the Ore completion of the set K[x\ \ {rUp} over shift and (/-shift algebras will 
contain zero. Thus, there is no analogon to the situation with (g)-Weyl algebras. 
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8.2. Open Problems 

1. From the description and the proof of the Algorithm 3.12, wc need in principle only 
the existence of a terminating Grobner basis algorithm over R resp. i?*. The latter relies 
on the constructivity of basic operations with fractions in R via i?*. Thus the Algorithm 
3.12 immediately extends to the case of a general OLGA {R,B,I) (cf. Dcf. 2.2) for a 
nontrivial completely prime ideal / C 5 by Bueso et al. (2003). In which bigger class 
of algebras one can perform concrete computations? One can work with local rings like 
^^Nmp, p G Jf" and K[[x]] for charii' = or with K{x} for K = R,C as coefficient 
domains. 

2. The algorithm computes U, V, D such that U MV = D. It is clear, that U, V are 
not unique. By fixing M and D, how can one describe the set {(i7, V) | UMV = D}1 
Does there exist U' resp. V from this set, such that special properties like unimodularity 
hold? This is connected, in particular, to the recent results of A. Quadrat and T. Cluseau 
(Cluzeau and Quadrat, 2010). Namely, they prove the existence of matrices U and V, 
which are unimodular over R^ for certain situations and present an algorithm for the 
computation of these matrices. 

3. One of the most important problems in non-commutative computer algebra is to 
show algorithmically, that two given finitely presented modules are not isomorphic. It is 
still open even for cyclic modules over general simple Euclidean Ore domain. Namely, 
let R be the latter domain and a,b & R with dega = deg6 > 0. Up to now there is no 
algorithm, which determines that R/{a) ^ R/{b). However, there are several situations, 
where this question has been solved. 

In Cluzeau and Quadrat (2008), the authors presented a semi- algorithm for finitely 
presented modules. Given degree bounds for the variables, the algorithm first looks for 
matrices, determining a homomorphism. If such a pair has been found, the kernel and 
cokernel of the corresponding homomorphism are computed and returned. 

4. As we have seen, over non-simple domains we have many different normal forms. 
Starting with classical operator algebras like the (g-)shift algebra and g-Weyl algebra, it 
is important to describe possible normal forms. 

5. As in the previous item, what is the most probable normal form for a random square 
matrix (which is then of full rank) over a non-simple domain? Though it seems that this 
might depend on the algebra, we conjecture the Diag(l, .., l,p), that is a Jacobson form, 
is the most probable one. Such approach might lead to the classification as in 4. 
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9. Appendix. The Code for an Example 

Consider the code for Example 3.15. After starting a Singular session, we need to 
set up the algebra first. 

> LIB "jacobson. lib" ; // load the library 

> ring w = 0, (x,s) , (a(0,l) ,Dp) ; 

The ring w is a commutative one in variables x,s over the field Q (0 stands for the 
characteristic). The substring (a(0,l) ,Dp) defines a monomial ordering on w, which in 

this case is an extra weight ordering, assigning weights to x and 1 to s. If two mono- 
mials are of the same weight, they are further compared with Dp (degree lexicographic 
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ordering). The described ordering mimicks the ordering on the rational shift algebra. For 
computations one can also use other orderings, like Dp or dp, cf. the Singular manual. 

> def W=nc_algebra(l , s) ; // set up shift algebra 

> setring W; 

This code creates and sets active a new ring W as a non-commutative ring from w with 
the relation s ■ x = 1 ■ xs + s. By executing W; a user will obtain the description of the 
active ring. 

> matrix m[2] [3] = 

> x*s-s+x"2-x, x*s+x"2, x*s+2*s+x"2+2*x, 

> s+x,0,s; 

> list J=jacobson(ni,0) ; 

Here the matrix was entered and the algorithm called. Note, that putting a string 
pr intlevel=2 ; before the j acobson call will output the progress of the algorithm. Higher 
values of printlevel will lead to more details printed during the execution. The list J 
has three entries, namely U, D, V. 

> print(J[l]*m*J[3]-J[2]) ; // check UinV=D 
==> 0,0,0, 

0,0,0 

> print (J [2]); // that is D 
==> 0,x4+3x3-x2-3x,0, 

0,0, X 

> print (J [1]); // that is U 
==> 0, -l/4x2-x-3/4, 

l/4,-l/4x-l/2 

> print (J [3]); // that is V 

==> ,_[1,2] ,xs2-s2+x2s-xs, 

_[2,1] ,_[2,2] ,_[2,3] , 

_[3,1] ,_[3,2] ,_[3,3] 
The symbols _ [2 , 2] before are printed, when the entries are long. One can access them 
as single polynomials 

> print (J [3] [2,2]) ; 

==> -21s2-7x2s-2xs+lls-3x3-2x2+13x-8 
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